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: ABSTRACT 

> 

CO , 

! We perform a systematic study of the Standard Model embedding in a D-brane config- 

CN . 

uration of type I string theory at the TeV scale. We end up with an attractive model and 

(Sj ; we study several phenomenological questions, such as gauge coupling unification, proton 

stability, fermion masses and neutrino oscillations. At the string scale, the gauge group is 

U(3) co i or x U(2) weak x U(l)i x U(l)b u ik. The corresponding gauge bosons are localized on 

^| three collections of branes; two of them describe the strong and weak interactions, while the 
r- 1 \ 

last abelian factor lives on a brane which is extended in two large extra dimensions with a 
size of a few microns. The hypercharge is a linear combination of the first three U(l)s. All 
remaining U(l)s get masses at the TeV scale due to anomalies, leaving the baryon and lepton 
numbers as (perturbatively) unbroken global symmetries at low energies. The conservation 
of baryon number assures proton stability, while lepton number symmetry guarantees light 
neutrino masses that involve a right-handed neutrino in the bulk. The model predicts the 
value of the weak angle which is compatible with the experiment when the string scale is 
in the TeV region. It also contains two Higgs doublets that provide tree-level masses to all 
fermions of the heaviest generation, with calculable Yukawa couplings; one obtains a natu- 
rally heavy top and the correct ratio mb/m T . We also study neutrino masses and mixings in 
relation to recent solar and atmospheric neutrino data. 



1 On leave of absence from CPHT, UMR du CNRS 7644, Ecole Polytechnique, 91128 Palaiseau, France 



1. Introduction 

In a previous work PJI2], a mm i ma l embedding of the Standard Model (SM) was proposed 
in a D-brane configuration of type I string theory with large internal dimensions and low 
fundamental scale 011]. The 577(3) color and SU(2) weak gauge fields were confined on 
two different collections of branes. The model correctly accommodated the right value of the 
weak angle for a choice of the string scale of a few TeV. It contained two Higgs doublets and 
guaranteed proton stability Among the issues, which were not addressed, are the fermion 
masses, neutrino oscillations, and a natural suppression of lepton number violating processes. 

A generic feature of the models studied was that some of the SM states should correspond 
to open strings with one end in the bulk, implying the existence of some extra branes, in 
addition to the ones used above [TJ Starting from the last point, in the present work 
we introduce an extra brane in the bulk with a corresponding U(l)b bulk gauge group 
0. This group is broken by anomalies, leaving behind an additional global symmetry that 
will be identified with the lepton number. In order to give masses to the neutrinos, we 
introduce a right-handed neutrino in the bulk j3] that carries non-trivial lepton number. 
Large neutrino masses are then forbidden by symmetry, while the right-neutrino coupling 
suppression required to explain the neutrino oscillation data, is achieved if the bulk has two 
dimensions of submillimeter size. 

More precisely, in the minimal case of one bulk neutrino, we show that solar and at- 
mospheric neutrino data can be accommodated using essentially the two lowest frequencies 
of the neutrino mass matrix: the mass of the zero mode, arising via the electroweak Higgs 
phenomenon, which is suppressed by the volume of the bulk, and the mass of the first 
Kaluza-Klein (KK) excitation. The former is used to reproduce the large mixing angle 
(LMA or even LOW) solution to the solar neutrino anomaly, through u e <-> transitions. 
The later is used to explain atmospheric neutrino oscillations with an amplitude which is 
enhanced due to logarithmic corrections of the two-dimensional bulk 0. Compatibility of 
the two conditions using one bulk right neutrino is possible only if one introduces a non- 
orthogonal angle between the two compact bulk dimensions, that leads simultaneously to a 
CP violation in the neutrino sector. Atmospheric oscillations contain however a significant 
sterile component which seems to be in contradiction with recent atmospheric data analyses. 

We also compute the tree-level Yukawa couplings of the two higgses to the fermions of the 
heaviest generation. They are given in terms of the gauge couplings and lead to a naturally 
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heavy top and a ratio rrib/m T compatible with the experimental data. Next, we proceed to 
a systematic description of the main features that we will use in the following sections. 

The general framework is type I string theory. We shall restrict ourselves to models in 
which the closed string sector is super symmetric, while supersymmetry is generically broken 
by the open strings at the string scale [7] . 2 Within our framework, the minimal ensemble of 
D-branes needed in our construction is the following mutually orthogonal stacks: a stack of 
three coincident branes to generate the color group, a second stack of two coincident branes 
to describe the weak SU{2) L gauge bosons, and one more brane to generate the U{l)b bulk 
discussed above. The resulting gauge group so far is U(3) c xU(2)lX Z7(l)&, with the three 
U(l) generators denoted by Q c , Ql and Qb, respectively. To ensure proton stability, we 
require baryon number conservation with generator B = Q c . The hypercharge Y cannot 
have a component along Q b , since this would lead to unrealistically small gauge coupling, 
and as explained in pQ the correct assignment of SM quantum numbers requires the presence 
of an extra abelian factor, named U(l)\ with generator Qi, living on an additional brane. 
This brane should lie on top of the color or the weak stack of branes, as we argue below. 

Since in our framework, supersymmetry is broken by combinations of (anti)branes and 

orientifolds which preserve different subsets of the bulk supesymmetries, any pair of D- 

branes Dp and Dp' satisfy p — p' = mod 4. It follows that a system with three stacks 

of mutually orthogonal branes in the six-dimensional internal (compact) space consists, up 

to T-dualities, of D9-branes with two different types of D5-branes, extended in different 

directions. Specifically, the £/(l)& lives on the D9-brane, while the U(3) c and U{2)l are 

confined on two stacks of 5-branes, the first along say the 012345 and the other along the 

012367 directions of ten-dimensional spacetime. Thus, the (submillimeter) bulk is necessarily 

two-dimensional (extended along the 89 directions), and the additional U(l)i brane has to 

coincide with either U(3) c or U(2)l. The parameters of the model are the string scale M s , 

the string coupling g s and the volumes v^, v§j and v$9 of the corresponding subspaces, in 

string units. 3 In terms of those, the four-dimensional Planck mass Mp is given by 

g 

M 2 P = —v 45 v G7 v 89 M* (1.1) 
9s 

and the non-abelian gauge couplings are 



11 11 . . 

= — W 4 5 ; ~2 = — v 67 (1-2) 

9t 9s 92 9s 



2 Recent progress in constructing type I vacua with structure close to the SM can be found in |51I§1ITU|. 
3 Using T-duality, we choose all internal volumes to be bigger than unity, > 1. 
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It follows that 

M 2 P = -^v 89 M* = —v m M 2 s , (1.3) 

gm «3«2 

where = gl/ATt and -Osg = f8g/(27r) 2 = R%R§ for a rectangular torus of radii R$,R9- The 
U(l)x gauge coupling gi is equal to (#2), if the t/(l)i brane is on top of the U(3) c (U(2)l). 

Upon T-duality, one finds two additional realizations: (i) a set of D3-branes (along 
0123) describing U(3) c , and two orthogonal sets of D7-branes along 01236789 and 01234567 
describing U(l)b and U(2) L , respectively; (ii) three sets of D5-branes along 012389, 012345 
and 012367, giving rise to U(l) b , U(3) c and U(2) L , respectively. In both cases, relation (jl.3|) 
remains intact. 

The gauge coupling g b of the U(l) b gauge boson which lives in the bulk is extremely 
small since it is suppressed by the volume of the bulk f 89 . For instance, in the case where 
the U(l) b lives on a D9-brane, its coupling is given by 

_ = _ W89 = __ (1 .4) 

where in the second equality we used eq. (jl.lj) . Using now the weak coupling condition 
g s < 1 and the inequality g s > gi 2 following from > 1 in eq. (|1.2|) , one finds 

/zM a V8M S 

v8— <^< — — , 1.5 

M P g 3 M P 

which implies that g b ~ 10" 16 - 10" 14 for M s ~ 1 - 10 TeV. 

If the U(l)b gauge boson is light, it will be subject to strong contraints coming from 
supernova observations, since it would be copiously produced in various nuclear reactions 
leading to supernova cooling through energy loss in the bulk of extra dimensions. The 
corresponding process is much stronger than the production of gravitons because of the 
non-derivative coupling of the gauge boson interaction ^T]. In fact, in the case of n large 
transverse dimensions of common radius R, satisfying m^,/? -1 «T with the gauge 
boson mass and T the supernova temperature, the production rate Pa is proportional to 

P A ~ gl x [R{T - m A )T x — ~ — , (1.6) 

where the factor [R(T — m^)] n counts the number of Kaluza-Klein (KK) excitations of the 
U(l)b gauge boson with mass less than T. This rate can be compared with the corresponding 
graviton production 



showing that for n = 2 (sub)millimeter extra dimensions, it is unacceptably large, unless the 
bulk gauge boson acquires a mass ^ 10 MeV. 

The paper is organized in seven sections, of which this introduction is the first. In 
Section 2, we perform a systematic search for models with four sets of branes corresponding 
to the gauge group U(3) c x U(2)l x U{1)i x U(l)b with the minimal standard model fermion 
spectrum and a Higgs sector that generates masses for all quarks and leptons of the heaviest 
generation. We identify the hypercharge U(l)y combination and in Section 3 we perform a 
renormalization group analysis of gauge couplings to identify models with low string scale, 
where the U(l)i is on top of either the color or the weak branes. In Section 4, we select four 
models with string scale in the TeV region, possessing in addition baryon and lepton number 
conservation, and we describe their main phenomenological features. 4 They all contain two 
Higgs doublets that can provide tree-level masses to all fermions. Moreover, apart from 
the hypercharge, all other abelian factors are broken by mixed gauge and gravitational 
anomalies and become massive at the string scale. In Section 5, we compute the tree-level 
Yukawa couplings of the two higgses to the fermions of the heaviest generation and study 
predictions for mass relations. In Section 6, we introduce one right-handed neutrino in the 
bulk and study the generation of neutrino masses and neutrino oscillations. Finally, Section 
7 contains our summary and conclusions. 

2. Model search 

As shown in pQ, the minimal D-brane configuration that can successfully accommodate 
the Standard Model (SM) consists of three sets of branes with gauge symmetry U(3) c x 
U(2) L x Z7(l)j_. The first set contains three coincident branes ("color" branes). An open 
string with one end attached to this set transforms as an SU(3) C triplet (or anti-triplet), but 
also carries an additional U(l) c quantum number which can be identified with the (gauged) 
baryon number. Similarly, U(2) L is realized by a set of two coincident branes ("weak" 
branes) and open strings attached to them from the one end are SU{2)i J doublets char- 
acterized by an additional U(1) L quantum number, the (gauged) weak "doublet" number. 
Moreover, consistency of the SM embedding requires the presence of an additional U(l) 1 
factor, generated by a single brane. This is needed for several reasons: TeV scale unification, 
baryon number conservation, and mass generation for all quarks and leptons of the heavi- 
4 Orientifold models with baryon and lepton number conservation were also constructed in Rcf. 9 . 
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est generation. The hypercharge is then a linear combination of the three abelian factors, 
y = ks Q c + &2 Ql + ki Qi, where Q c , Ql,Qi are the charges under U (1) c , U (1) L , U (1) 1 re- 
spectively. It turns out pQ that there exist four possible "viable" models that reproduce the 
weak mixing angle al low energies. They correspond to k^ — | (k 3 = — ~), /c 2 = ±|, ki = 1 
and require the abelian brane 11(1)-^ to be on top of the color (weak) branes, so that g 3 = g± 
(92 = 9x)- 

In all the above brane configurations there exist states (e.g. the SU(2)i singlet anti- 
quarks) which correspond to open strings with only one of their ends attached to one of the 
three sets of D-branes. The other end is in the bulk, and requires the existence of some 
additional branes extended in the bulk, carrying extra quantum numbers. In this work, we 
consider a minimal extension of the models considered in pQj by introducing one additional 
D-brane in the bulk giving rise to an extra abelian gauge factor U(l) b . As we will see later, 
the requirement of baryon and lepton number conservation leads to four possible models that 
we are going to study in the next section. However, in this section, we do not impose this 
constraint and we systematically explore the possibility of reproducing the SM spectrum, 
together with possibly additional Higgs scalars, as open strings stretched between any two 
of the four sets of branes. The extension of the Higgs sector is required for the realization 
of the electroweak symmetry breaking and mass generation for all fermions of at least one 
(the heaviest) generation. 

Thus, the total gauge group is 

G = U(3) c xU(2) L xU(l) 1 xU(l) b 

= SU(3) c xU(l) c xSU(2) L xU(l) L xU(l) 1 xU(l) b (2.1) 

and contains four abelian factors. The assignment of the SM particles is partially fixed from 
its non-abelian structure. The quark doublet Q corresponds to an open string with one end 
on the color and the other on the weak set of branes. The anti-quarks u c , d c must have one 
of their ends attached to the color branes. The lepton doublet and possible Higgs doublets 
must have one end on the weak branes. However, there is a freedom related to the abelian 
structure, since the hypercharge can arise as a linear combination of all four abelian factors. 
In a generic model, the abelian charges can be expressed without loss of generality in terms 
of ten parameters displayed in Table [TJ 

In a convenient parametrization, normalizing the U(N) ~ SU(N) x U(l) generators 
as TrT a T b = 5 ab /2, and measuring the corresponding U(l) charges with respect to the 
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particle 


U(l) c 


U(1)l 


U{1\ 




Q(3,2, |) 


+1 


w 








u c (3,l,-D 


-1 







a 2 


d c (3,l,+|) 


-1 





bi 


b 2 


MM,-!) 





+1 


Cl 


c 2 


e c (l,l,+l) 





d L 




d 2 



Table 1: SM particles with their generic charges under the abelian part of the gauge group 
U(3) c x U(2) L x 1/(1)! x U{1\. 

coupling g/y2N, the ten parameters are integers: ai 5 2, &i,2, ci,2, d 2 = 0, ±1, d\ = 0, ±1, ±2, 
<7l = 0, ±2, w = ±1 satisfying 

En = En = En = 1 ' E i*i = 2 - ( 2 - 2 ) 

i=l,2 i=l,2 1=1,2 i=l,2,i 

The first three constraints in ()2.2)1 correspond to the requirement that the u c and d c anti- 
quarks, as well as the lepton doublet, must come from open strings with one end attached 
to one of the abelian D-brane sets. The fourth constraint forces the positron e c open string 
to be stretched either between the two abelian branes, or to have both ends attached to the 
abelian U(l) l brane, or to the weak set of branes. In the latter case, it has U(1)l charge 
±2 and is an SU(2) l singlet arising from the antisymmetric product of two doublets. The 
parameter w in Table^refers to the U(1) L charges of the quark- doublets, that we can choose 
to be ±1, since doublets are equivalent with anti-doublets. Note that a priori one might also 
consider the case in which one of the u c and d c anti-quarks arises as a string with both ends 
on the color branes (3 x 3 = 3 + 6), so that its U{l) c charges would be ±2. This, however, 
would invalidate the identification of U(l) c with the baryon number and forbid the presence 
of quark mass terms, since one of the combinations Qu c and Qd c would not be neutral under 
U{l) c . Hence, this case will not be explored. 

The hypercharge can in general be a linear combination of all four abelian group factors. 
However, we restrict ourselves to models in which the bulk U(l) b does not contribute to the 
hypercharge, in order to avoid an unrealistically small gauge coupling. Hence, 

Y = k 3 Q c + k 2 Q L + k 1 Q 1 . (2.3) 
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The correct assignments for SM particles are reproduced, provided 

k 3 + k 2 w = - 
o 

2 

-k 3 + a t k x = -- 

-h + b 1 k 1 = ~ (2.4) 

k 2 + ci ki = -- 
k 2 d,L + d\ki = 1. 

Notice that the second and third of the above equations imply that k\ ^ 0. 

The next step, after assigning the correct hypercharge to the SM particles, is to check 
for the existence of candidate fermion mass terms. Here, we discuss only the question of 
masses for one generation (the heaviest) and we do not address the general problem of 
flavor. To lowest order, the mass terms are of the form Qd c H\ L , Qu c H u and Le c H\ where 
H d , H u , H e are scalar Higgs doublets with appropriate charges. In a generic model, there 
are four different candidate Higgs scalar doublets (and their conjugates) Hi, . . . ,H±, with 
U(1) L x U{1) 1 x U{l) b charges: 

{Hi, H 2 , H 3 , H,} = {(1, 1, 0), (1, 0, 1), (1, -1, 0), (1, 0, -1)} . (2.5) 

It is easy to show that for any hypercharge embedding of the form ()2.3|) with ki ^ 0, there are 
at most three of the above Higgs doublets that have the correct hypercharge. Depending on 
the parameters of the model, they can be reduced to two. For the generic charge assignments 
of Table [TJ the required Higgs charges are 

H u = (1,2,0, -w,-ai,-a 2 ) 

H d = (1, 2, 0, +w, +h, +b 2 ) (2.6) 
H e = (1, 2, 0, 1 + d L , ci + di, c 2 + d 2 ) . 

Provided the constraints (J2.2|) are satisfied, both H u and H d have the right charges of 
(|2.5|) and correspond to strings stretched between the weak and one of the abelian branes. 
Thus, ()2.2|) guarantees the existence of tree-level quark masses. On the other hand, the 
existence of H e depends on the particular choice of parameters, e.g. for ci + di = 2, H e does 
not exist and a tree-level lepton mass term (Le c H^) is forbidden. The generic constraint 
that guarantees tree-level lepton masses is 

\ci + dt\ = \l + d L \ = 1 . (2.7) 

i=l,2 
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Given the smallness of the lepton mass compared to the masses of the quarks, of the same 
generation, it would be reasonable to examine also the possibility that the lepton mass is 
generated by a higher order term. The next order candidate lepton mass term is of dimension 
six, proportional to Le c WH'H. The constraint in this case is more complicated and the 
method we are going to use is the following: for each configuration that satisfies all other 
constraints except ([2.7)1 . we derive explicitly the candidate Higgs doublets and check the 
existence of possible fifth order mass terms. 5 

The hypercharge constraints (|2.4|) can be easily solved. They require a\ 7^ hi and 

h = itjt 7 (2-8) 

3 (61 - a x ) 

k, = - 'f^ (2.9) 

2 (Oi — <2i) w 

h = (2.10) 
b\ — a\ 

ci = _^i + K±M (2 . n) 

1 2 2w K J 

a h , ( Q i + b i) d L 

di = bi-ai-\ . (2.12) 

2w 

The allowed values of (a x , 61) are {(-1, 0), (-1, 1), (0, -1), (0, 1), (1, -1), (1, 0)}. However, 
we notice that the solutions with parameters (ai, 61, ci, di, ki) and (— Oi, —61, — Ci, —di, —k\) 
are equivalent, since they correspond to a global change of sign Qi — ► —Qi. Thus, it is 
sufficient to search for solutions with (a 1; foi) G {(— 1, 0), (— 1, 1), (0, 1)}. Solving for these 
choices, we get three allowed hypercharge embeddings: 

(i) o x = -l,6i = l : Y = \Q c + l -Qi (2.13) 

o 2 

I w 

{11) 0l = -1,61 = : r = --g c + -g L + gx (2.14) 

2 10 

(m) ai = 0,6i = l : Y = -Q c --Q L + Qi. (2.15) 

Case (i) leads to C\ = — l,c 2 = 0,di = 2,6?! = d L = 0. This is a special solution where the 

U(l) b brane decouples from the model since no SM particles are attached to it. It satisfies 

(|2.7j) and thus leads to tree level lepton masses. The solution exists for both w — ±1, 

as the value of w does not play an important role when k 2 = 0. In case (ii), we have 

Here, we check only the conservation of all gauge quantum numbers. In the string context, there may 
be additional selection rules for the non- vanishing of the corresponding couplings that are model dependent. 
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c\ = — (1 +w)/2, d,L = 0, d± = 1 or c\ = (1 +w)/2,cIl = 2w,di = d 2 = 0, while case (iii) 
leads to c\ = (w — l)/2, dx = 0, cfi = 1 or c\ — (1 + w)/2, di = 2w, d\ = d 2 = 0. 

Combining the above three cases with the constraints ()2.2j) and ()2.7j) . we get 9 distinct 
configurations with tree-level quark and lepton masses, displayed in the upper part of Ta- 
ble El Relaxing the constraint ()2.7|) with the requirement that lepton masses arise through 
dimension six effective operators, leads to 6 more distinct models corresponding to the cases 
10-15 of Table 121 In deriving these configurations, we have eliminated all models connected 
to the ones above by the global charge redefinition Q b — > —Qb- 

As we mentioned before, in all the above configurations, we can define the baryon number 
B as 

B = ^Q C . (2.16) 

As we will argue below, U(l) c gauge invariance is broken by anomalies to a global symmetry, 
implying baryon number conservation in type I string perturbation theory Since lepton 
number is also conserved at present energies, we can further examine which of the above 
models possess also the lepton number L as a symmetry. In general, L can also be expressed 
as a linear combination of all abelian factors, 



ViQi (2.17) 



=c,L,l,i 



that satisfies 

Pc + Plw = 

-p c + aipi + a 2 p b = 

-p c + bxpi + b 2 p b = (2.18) 

PL + Cipx + C 2 Pb = 1 

d L PL + d l p l + d 2 p b = -1 

Inspection of (|2.18|h in conjunction with ()2.4j) that requires a± ^ b\, implies that lepton 
number can only be defined for p b ^ 0, i.e. only in the presence of the bulk U(l)b- This 
is of course expected, since the models without U(l)b have no lepton number |l r Solving 
explicitly ()2.18|) for each one of the cases of Table [21 we find that only four models, namely 
2,4,6,9, incorporate the lepton number as a (gauged) abelian symmetry. Its precise definition 
for each of these models is also presented in the table. 
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1 

2 
3 
4 
5 
6 
7 
8 
9 



ai 
-1 
-1 
-1 





-1 
-1 







a 2 




1 
1 




1 
1 



1 




1 
1 




1 
1 





-1 
-1 




-1 
1 






Cl 

-1 



-1 



-1 








c 2 



di 
2 
1 
1 
1 
1 







d 2 



1 

1 

1 

1 











d L 






-2 
-2 
-2 

_9 



W 

1 
-1 

1 

1 
-1 
-1 
-1 
1 
1 



Y 

\Qc + \ Qi 
—\ Qc — \ Ql + Qi 

-\Qc+\Q L + Ql 

I Qc — \ Ql + Qi 

I Qc + h Ql + Qi 



\Qc~ \ Ql + Qi 
\ Qc — \ Ql + Qi 



\ Qc — \ Ql + Qi 
\Q c -\Ql + Qi 



\ Qc + \ Ql — \Qi — \Qb 



\Qc + \ Ql — \Qi — \ Qb 



I Qc + k Ql — \Qi — \Qb 



■\Q c + \Q L -\Qi-\Qb 



10 -1 1 -1 2 1 -\Q C + \Q L + Qi 

11 1 1 -1 2 -1 IQc + \Ql + Qi 

12 -10 1 1 1 1 -1 -\Q c -\Ql + Qi 

13 -1 -1 1 1 1 -1 -\Q c -\Ql + Qi 

14 1 1 1 1 1 1 IQ c -\Ql + Qi 

15 o li o o-ii-io l \Q c -hQ L + Qi 



Table 2: Distinct models with lepton masses generated either at tree level Le c H^ (cases 1-9), or by dimension six effective 
operators ~ Le c (W H) / 'M 2 (cases 10-15). We also display the lepton number combination L (when it exists) and the 
number of Higgs doublets nh, needed to generate quark and lepton masses. 



3. The weak angle and the string scale 



We now come to the determination of the string scale consistent with the low energy 
SM data. Following the hypercharge definition (|2.3|) . the low energy data depend on the 
couplings g 3 , g 2 and gi of the three brane sets U(3) c , U{2) L and U(l) v These couplings are 
in principle independent, but, as already explained in the introduction, in order to lower the 
string scale we have to consider configurations where the U(l) 1 brane is on top of either the 
U(3) c or the U(2) L stacks. Hence, we have two possible coupling relations at the string scale 

(*) 93 = 9i or (ii) g 2 = gi ■ (3.1) 

In our normalizations, the hypercharge coupling g Y at the string scale is expressed as 

1 6fcf 4Jfc| 2k\ 
9y 95 92 91 

Following the one loop coupling evolution (ctj = gf/An), 

1 1 bi A*M S 

^(Ms) ~^{Mz) + 47 n ll^' (3 ' 3) 

where 63 = —7, b 2 = —10/3 + n^/6, by = 20/3 + n^/6 and nu is the number of scalar Higgs 
doublets. The constant A 1 corresponds to a model independent piece of the type I string 
thresholds, entering in the identification of the string scale with the ultraviolet cutoff of the 
effective field theory. Its value was computed in Ref. ^2j to be A 1 = 1/ ' ~ 0.4, where 7 
is the Euler's constant. When the string scale is very high compared to present energies, this 
represents a small correction compared to the dominant logarithmic contribution coming 
from the renormalization group evolution. On the other hand, when the string scale is low, 
such a correction becomes important, as it effectively changes the string scale by roughly 
a factor of two, and should be taken with caution since it is of the same order with the 
(unknown) model dependent part of threshold corrections. Consequently, we will leave A 1 
as a parameter and discuss its possible effects on our results case by case. 

Solving the one-loop renormalization group equations (RGE) for the coupling evolution, 
the values of the weak mixing angle sin 2 9w and of the strong coupling 03 at the Z-mass Mz 
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are related to the couplings at the string scale: 

sin W (M Z ) = : I 7- r^ ln ^rr ^ (3-4) 

v 7 l + k Y 2?r (l + M M z v ' 



a em (M z ) 
1 + k 



Y 



6k l ( — 7Y7\ 7TFT J + 2A; i ' 



a L (M s ) « 3 (M S )7 1 V«l(M s ) ai (M 



a 3 (M z ) " a em (M z ) 1 + A; y 2vr 1 + A; y ° g M z 1 ' j 



1 



1 + Jfe 



^a L (M s ) a 3 (M s ); \a L {M s ) a x {M t 

where ky = 6 fcf + 4 fc| + 2/c 2 and a em is the electromagnetic coupling. 

Given a coupling relation of ()3.1j) . we can use eqs. (J3.4)) and ()3.5j) to determine the string 
scale M s that correctly reproduces the low energy data. Clearly, the solution depends on 
|/c 3 |, \k 2 \ and \k\\. According to our previous analysis, there are three classes of models, 
which correspond to the three possible hypercharge embeddings (j2.13|) . ()2.14j) and ([2.15)1 : 
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Using ()3.4j) and (|3.5jl . for each of the embeddings ()3.6|) and the unification conditions (|3.1|) . 
we computed the corresponding string "unification" scale My = A 1 M s . In our calculation we 
have used the following values for the low energy quantities a 3 (M z ) = 0.119, sin 2 6w = 0.231, 
a em (M z ) = 1/127.934. The results are presented in Table El 

In the above calculations we have assumed that the number of doublets nh is the minimum 
rih = 2 required by the model. Of course, one can consider models with more doublets which 
can be for instance replicas of these two. It would be thus interesting to examine the 
dependence of the above results on the number of doublets. To this end we can extract 
analytic formulas regarding the unification scale My. For the case g± = g 3 , taking for 
simplicity k% — 1 and ky = ±|, we find 

3(4 + 7fci) log ^ = -4tft(1 - 2 sin 2 MM*)) " 2(1 + 3*2)-^ (3.7) 
7T M z a em {M z ) a 3 {M z ) 

which implies that at the one-loop My is independent of the number of doublets. Similarly, 

for #1 = g 2 and k\ = 1, k L = ±|, we have 
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Table 3: The string unification scale My and the two independent gauge couplings for the 
two possible brane configurations and the various hypercharge embeddings. 

where we find a very weak dependence. Obviously, the number of doublets affects the value 
of the weak gauge coupling at M s and thus the volume of the bulk through (jl.3j) . 

4. The models 

So far, we have classified all possible Z7(3) c x U(2) L x U(l) 1 x U(l) b brane models that 
can successfully accommodate the SM spectrum. The quantum numbers of each model as 
well as the hypercharge embedding are summarized in Tabled Furthermore, compatibility 
with type I string theory with string scale in the TeV region, requires the bulk to be two- 
dimensional of (sub) millimeter size, and leads to two possible configurations: Place the U(l) 1 
brane on top of the weak U(2) L stack of branes or on top of the color Z7(3) branes. These 
impose two different brane coupling relations at the string (unification) scale: g\ = g 2 or 
gi = g 3 , respectively. For every model, using the hypercharge embedding of Table the 
one loop gauge coupling evolution and one of the above brane coupling conditions, we can 
determine the unification (string) scale that reproduces the weak angle at low energies. The 
results are summarized in Table El 

According to the results of Section 2, there are three distinct hypercharge embeddings 
that correspond to (|Jfe 3 |, |Jfe 2 |, N) = {(1/6, 0, 1/2), (1/3, 1/2, 1), (2/3, 1/2, 1)}. Since we wish 
to restrict ourselves to models in which supersymmetry is broken at the string scale (M s ), 
we would like M s to be low, at the TeV scale, to protect the mass hierarchy. Thus, model 
1 of Table El with hypercharge embedding (1/6,0,1/2), is rejected for both U(l) 1 brane 
arrangements, since the resulting string scale is too high. Furthermore, models with hyper- 
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charge embedding (1/3, 1/2, 1) lead to M s ~ 10 3 TeV for g\ = g%. This scale, although much 
lower than the traditional GUT scale, is rather high for the stabilization of hierarchy. On 
the contrary, for g\ = g2, we find M s ~ 0(1) TeV (when the universal threshold correction 
A 1 is taken into account) that lies at the edge of the present experimental limits. The third 
embedding (2/3, 1/2, 1) reproduces successfully the low energy data only for gi = g 3 and a 
string scale M s ~ O(10) TeV. 

In all configurations of Tabled the baryon number appears as a gauged abelian symmetry. 
This symmetry is broken due to mixed gauge and gravitational anomalies leaving behind 
a global symmetry. Baryon number conservation is essential for low string scale models, 
since one needs to eliminate effective operators to very high accuracy in order to avoid 
fast proton decay, starting with dimension six operators of the form QQQL which are not 
sufficiently suppressed fT5] . In addition to baryon number, one should also assure that the 
lepton number is a good symmetry of the low energy theory. Lepton number conservation 
is also essential for preservation of acceptable neutrino masses, as it forbids for instance the 
presence of the dimension 5 operator LLHH. Such an operator would lead to large Majorana 
neutrino masses, of the order of a few GeV, in models where the string scale, typically a 
few TeV, is too low for the operation of an effective sea-saw mechanism. Hence, we shall be 
interested only in models in which the lepton number is a good symmetry. Indeed, as seen 
in Table EJ only in four models, namely 2,4,6 and 9, lepton number appears as a (gauged) 
abelian symmetry. Being anomalous, this symmetry will be broken, but lepton number will 
survive as a global symmetry of the effective theory. 

In fact, these four models can be derived in a straightforward way by simple considerations 
of the quantum numbers. The quark doublet Q is fixed by non abelian gauge symmetries, 
while existence of baryon number implies that the antiquarks u c , d c correspond to strings 
stretched between the color branes and one each of the abelian branes U{1)\ and £7(1)&. 
Thus, one has two possibilities leading to models that we call A (d c has one end in the bulk) 
and B (u c sees the bulk). Existence of lepton number fixes the lepton doublet as a string 
stretched between the weak branes and the U(l)b brane, while for each of the models A and 
B there are two possibilities for the antilepton e c to emerge as a string stretched between the 
two abelian branes, or to have both ends on the weak branes. Thus, we obtain two additional 
models that we call A' and B'. As it can also be seen in the table, all these models have 
tree-level quark and lepton masses and make use of only two Higgs doublets. They also 
require low energy string scale for some of the brane coupling conditions. We now proceed 
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to a detailed study of these four models and to an analysis of their main phenomenological 
characteristics. 

Notice from Table El that in both classes of models A and B, the coupling constant ratio 
is gilds — 0-6 at the string scale, implying through the relations of Section 1 that at least 
one of the internal compact dimensions along the world-volume of the weak set of branes 
must be larger than the string length, by at least a factor of two (in the case of two large 
dimensions, or by a factor of four in the case of one). The relevant experimental signal would 
be the production of Kaluza-Klein excitations for the W ± bosons and the other mediators 
of the electroweak interactions but not of gluons, providing one of the first indications of 
new physics |T3] . 

Models A and A' 

We consider here the models 2 and 6 of Table |2l hereafter referred as models A and A' 
respectively. They are characterized by the common hypercharge embedding 

y = -\qc-\q l + Qi (4.i) 

but they differ slightly in their spectra. The spectrum of model A is 

Q (3, 2, +1,-1, 0,0) 
u c (3, 1,-1, 0,-1,0) 
cf(3, 1,-1, 0,0,-1) 
L(l, 2,0, +1,0,-1) 
e c (l,l,0,0,+l,+l) 
#„(1,2,0,+1,+1,0) 
^(1,2,0,-1,0,-1) 

while in model A' the right-handed electron e c is replaced by an open string with both ends 
on the weak brane stack, and thus e c = (1,1,0,-2,0,0). The two models are presented 
pictorially in Figure ^ 

Apart from the hypercharge combination (|4.1|) all remaining abelian factors are anoma- 
lous. Indeed, for every abelian generator Qi, I = (c, L, 1, b), we can calculate the mixed gauge 
anomaly Kjj = TtQjTj with J = SU(3), SU (2), Y, and gravitational anomaly Kj^ = Tr Qj 
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Figure 1: Pictorial representation of models A, A'. 
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(4.2) 



It is easy to check that the matrices KK T for both models have only one zero eigenvalue cor- 
responding to the hypercharge combination (|4.1j) and three non vanishing ones corresponding 
to the orthogonal U(l) anomalous combinations. In the context of type I string theory, these 
anomalies are cancelled by a generalized Green-Schwarz mechanism which makes use of three 
axions that are shifted under the corresponding U(l) anomalous gauge transformations [To] . 
As a result, the three extra gauge bosons become massive, leaving behind the correspond- 
ing global symmetries unbroken in perturbation theory [Jj]]. The three extra U(l)'s can be 
expressed in terms of known SM symmetries: 



Baryon number B 
Lepton number L 
Peccei-Quinn Q P q 



- {Qc + Ql-Qi 



(4.3) 



(Q c -Ql-3Q!-3Q 6 



Thus, our effective SM inherits baryon and lepton number as well as Peccei-Quinn (PQ) 
global symmetries from the anomaly cancellation mechanism. Note however that PQ is 
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the original Peccei-Quinn symmetry only in model A', such that all fermions have charges 
+1, while H u and Hd have charges —2 and +2, respectively. In model A, the global PQ 
symmetry defined in Q4.3J) is similar but with lepton charge +3. The reason is that in model 
A the fermion-Higgs Yukawa couplings are different, and leptons get masses from H u and 
not from Hd- 

The general one-loop string computation of the masses of anomalous U(l) gauge bosons, 
as well as their localization properties in the internal compactified space, was performed 
recently for generic orientifold vacua |17j . It was shown that orbifold sectors preserving 
N = 1 supersymmetry yield four-dimensional (4d) contributions, localized in the whole six- 
dimensional (6d) internal space, while N = 2 supersymmetric sectors give 6d contributions 
localized only in four internal dimensions. The later are related to 6d anomalies. Thus, even 
U(l)s which are apparently anomaly free may acquire non-zero masses at the one- loop level, 
as a consequence of 6d anomalies fTJ E| . These results have the following implications in 
our case: 

1. The two U(l) combinations, orthogonal to the hypercharge and localized on the strong 
and weak D-brane sets, acquire in general masses of the order of the string scale from 
contributions of N = 1 sectors, in agreement with effective field theory expectations 
based on 4d anomalies. 

2. Such contributions are not sufficient though to make heavy the third U(l) propagating 
in the bulk, since the resulting mass terms are localized and suppressed by the volume 
of the bulk. In order to give string scale mass, one needs instead N = 2 contributions 
associated to 6d anomalies along the two large bulk directions. In our models such 
contributions are indeed in general present and arise from mixed 6d gauge-gravitational 
anomalies of two different sources: (i) the generic presence of a neutral 6d Weyl fermion 
on the bulk brane which coincides either with the £/(l)& gaugino (in the supersymmetric 
case) or the goldstino of the non-linearly realised supersymmetry (in the brane SUSY 
breaking case |7j); (ii) the contribution of the right-handed neutrino which arises from 
a six- dimensional Weyl spinor. As a result, the third abelian gauge field U(l)b acquires 
also a mass of the order of the string scale, although its gauge coupling is tiny due to 
the volume suppression (see eq. (|1.5j) ). 

3. Special care is needed to guarantee that the hypercharge remains massless despite the 
fact that it is anomaly free, along the lines of Ref . ^7] . 
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The presence of massive gauge bosons associated to anomalous abelian gauge symmetries 
is generic. Their mass is given by M\ ~ g s M^, up to a numerical model dependent factor 
and is somewhat smaller that the string scale. When the latter is low, they can affect low 
energy measurable data, such as g — 2 for leptons jTU] and the p-parameter (201 , leading to 
additional bounds on the string scale. 

Note that the global PQ symmetry leftover from U(l)b is spontaneously broken by the 
Higgs expectation value giving rise to an unwanted electroweak axion. A possible way out 
was suggested in ref. [IJ, using an appropriate departure away from the orientifold point. 

A plausible extension of the model is the introduction of a right-handed neutrino in the 
bulk. A natural candidate state would be an open string ending on the U(l) b brane. Its 
charge is then fixed to +2 by the requirement of existence of the single possible neutrino 
mass term L v R . The suppression of the brane-bulk couplings due to the wave function of 
Ur would thus provide a natural explanation for the smallness of neutrino masses. Note that 
if the zero mode of this bulk neutrino state is chiral, the anomaly structure of the model 
changes: B — L becomes anomaly free and as a consequence the associated gauge boson 
remains in principle massless. However, as we discussed above, this is not in general true 
because of 6d anomalies [T7j. In any case, this problem is absent if we introduce a vector-like 
bulk neutrino pair 

Ml, 1,0, 0,0, +2) + 4(1, 1,0, 0,0, -2) 

that leaves the anomalies (|4.2|) intact. Note that u R does not play any role in the subsequent 
discussion of neutrino masses and oscillations. 

Coming to the issue of gauge couplings and the string scale, as already explained we have 
two different realizations for each model. The first is with g\ = g^ at M s that corresponds 
to a configuration where the Z7"(l) 1 brane is placed on top of the color branes. According 
to Table E3 this leads to an intermediate string scale M s ~ 10 6 GeV, which appears too 
high to guarantee the stabilization of hierarchy. The second possibility is to take the U(l) 1 
brane on top of the weak branes, leading to g± = #2- The required string scale is now 
low M s ~ O(500) GeV (300-800 GeV, depending on the threshold corrections), and could 
account for the stability of the hierarchy. 

Models B and B' 

Another phenomenologically promising pair of models consists of solutions 4 and 9 of 
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Figure 2: Pictorial representation of models B and B' . 

Table |2l named hereafter B and B', which corresponds to the hypercharge embedding 

Y = \Qc~\Ql + Qi. (4.4) 

The spectrum is 



Q(3,2,+l,+l,0,0) 
u c (3, 1,-1, 0,0,1) 
<f(3, 1,-1, 0,1,0) 
L(l, 2,0, +1,0,-1) 
e c (l,l,0,0,+l,+l) 
# u (l,2,0,-l,0,-l) 
tf,(l,2,0,+l,+l,0) 

for model B, while in B' e c is replaced by e c (l, 1, 0, —2, 0, 0). The two models are represented 
pictorially in Figure |21 The four abelian gauge factors are anomalous. Proceeding as in the 
analysis ()4.2}) of models A and A' , the mixed gauge and gravitational anomalies are 

/ 1 I I \ ( 



K (B) 



o i I §\ 

§20-1 

_3 2 4 11 
2 3 3 6 



\ 8 4 2 / 



1 i i\ 

x 2 2 







4 15 
'3 3 6 



\ 6 



(4.5) 



1 / 



20 



It is easy to see that the only anomaly free combination is the hypercharge (|4.4|) which 
survives at low energies. All other abelian gauge factors are anomalous and will be broken 
by the generalized Green-Schwarz anomaly cancellation mechanism, leaving behind global 
symmetries. They can be expressed in terms of the usual SM global symmetries as the 
following U(l) combinations: 

Baryon number B = — Q c (4.6) 

3 

Lepton number L = (Q c - Q L + Qi + Q b ) (4.7) 

Peccei-Quinn Q PQ = ^(-Q c + 3Ql + Qi + Qb) (4.8) 

Similarly to the analysis of models A and A', the PQ charges defined above are the traditional 
ones only for model B. In model B', the lepton charge result of the Higgs Yukawa 

couplings to the fermions (see below). The right handed neutrino can also be accommodated 
as an open string with both ends on the bulk abelian brane: 

^(1,1, 0,0,0, +2) + 4(1, 1,0, 0,0, -2) 

According to the RGE running results of Table El there is only one brane configuration, 
for the models under discussion, that reproduces the weak mixing angle at low energies. This 
consists of placing the U(l) 1 brane on the top of the color branes, so that g\ = g%, which 
leads to M s ~ O(10) TeV (7-17 TeV, depending on the threshold corrections). 

5. Fermion masses 

Although the general question of quark and lepton masses goes beyond the scope of this 
paper, we would like to make here some comments in the context of our constructions. The 
Yukawa couplings relevant to fermion masses are constrained by the various U(l) symmetries 
and can present interesting patterns. 

• Model A. The relevant Yukawa couplings are 

M A = \ u Qu c H u + \ d Qd c Hl + \ e Le c Hl + \„LH d u R (5.1) 

Here, charged leptons and up quarks (of the heaviest generation) obtain masses from the 
same Higgs (H u ). 
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When all Yukawa couplings arise at the lowest (disk) order, it is easy to check that in the 
simplest case (absence of discrete selection rules, etc), they satisfy the following relations: 

\ u = \ e = ^/2g 2 , X d = \/2g~s , K = V2g b . (5.2) 

The top and bottom quark masses are given by: 

m = g 2 v sin (3 ; m b = ^fgv cos /3 , (5.3) 

where tan (3 = v u / Vj, with v u and the vacuum expectation values ( VEVs) of the two 
higgses H u and Hd, respectively, and v = \Jv \ + v\ = 246 GeV. Note that in the case where 
the color branes are identified with D3 branes, one has ^fgl = g$, and in any case g s > g\. 
Note also that since the string scale in this model is relatively low, M s ^ 1 TeV, there is no 
much evolution of the low energy couplings from the electroweak to the string scale. Thus, 
using the known value of the bottom mass ~ 4 GeV, one obtains for the top quark mass 
m t ~ 162 GeV which is less than 5% below its experimental value mj Xp = 174.3 ±5.1 GeV. 
In addition, the Higgs VEV ratio turns out to be large, tan/3 ~ 100. Note that such a 
large value is not in principle problematic as in the supersymmetric case, but it can lead to 
important higher order corrections. 

On the other hand the r-mass is of the same order as the top mass, which is unrealis- 
tic. However, there is still the possibility that the lepton Yukawa coupling A e vanishes to 
lowest order due to additional string discrete selection rules, and is generated by a higher 
dimensional operator of the form Le c {H\ji^ H) providing the appropriate suppression. 6 

• Model A'. The Yukawa couplings here are 

M A , = \ u Qu c H u + \dQd c H d + \ e Le c H d + \ u LH d is R (5.4) 

with the same relation for the tree-level couplings as in (J5.2|) . Using the parametrization 

in ()5.3|) we see that the relation of m t to m& is the same as in model A and the same 

remarks apply. Since here the lepton and down quark acquire their masses from the same 

Higgs, one obtains the phenomenologically interesting relation: m h /m T = J~g~ s lg% = Qzl Qi, 

when strong interactions are on D3 branes. Thus, from Table El mb/m T ~ 1.75 at the 

(string) unification scale, which is in the upper edge of the experimentally allowed region 

6 Models with similar properties have been considered in the past in the perturbative heterotic string 
framework. 
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at the Z-mass, 1.46 ^ m b /m T \ exp ^ 1.75. This relation could replace the successful GUT 
prediction m& = m T of the conventional unification framework, in low scale string models. 
In conclusion model A' seems to be able to generate the required hierarchy of masses for the 
third generation. 

• Model B. The relevant trilinear Yukawa couplings are, 

M B = \ u Qu c H u + \ d Qd c Hl + \ e Le c Hl + \ v LH u u R (5.5) 
The tree-level Yukawa couplings satisfy 

\ e = \ u = y^2g~ s , X d = V2g 3 , X u = V2g b (5.6) 

and we have 

m t = y/g s vsmP ; m b = g 3 v cos (3 . (5.7) 

The first relation implies again a heavy top, while the bottom to tau mass ratio is now 
predicted, with a value m b /m T = g 3 j ^fgl ^ 1 which is apparently far from its experimental 
value. However, in this case, the string scale is relatively high and therefore one should 
take into account the renormalization group evolution above the weak scale. Solving the 
associated RGEs with the boundary conditions f)5.7j) and assuming #3 = y/gl, we obtain 
acceptable m b and m T masses for M s ~ 3 x 10 3 TeV and tan/? ~ 80. Note that the 
successful prediction of m b and m T is related to the condition m b = m T at the (string) 
unification scale, which in the case of non-supersymmetric Standard Model is obtained at 
relatively low energies |2T]. Indeed, in Figure 01 we plot the mass ratio m b /m T as a function 
of the energy, within the non supersymmetric Standard Model with two Higgs doublets. 
Nevertheless, the resulting value of M s is still significantly higher than the unification scale 
required from the analysis of gauge couplings in section 3. Moreover, the top quark mass 
turns out to be rather high, m t ~ 220 GeV. It is an open question wether this discrepancy can 
be attributed to threshold corrections that can be important in the case of two dimensional 
bulk jEj. 

• Model B'. The relevant Higgs couplings are given by 

M B , = X u Qu c H u + X d Qd c H\ + X e Le c Hl + X v LH u v R (5.8) 
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Figure 3: Evolution of the ratio m^jm T as a function of the energy fi for tan/? = 2 and 
tan (3 — 80. We have used as low energy parameters nib = 4 GeV, m top = 174 GeV a 3 (M z ) = 
0.12, sm 2 9 w = 0.23113. 

while the tree-level Yukawa couplings by 

X u = ^/2g~ s , X d = \/2g 3 , X u = \flg h and A e = \/2g 2 (5.9) 

Here, as in model A, the r and top mass are of the same order and thus in conflict with 
experiment. As in model A, vanishing leading order coupling could be a way out. 

In the above analysis we have also assumed that only the heaviest generation acquires 
masses at the lowest order. The other two are considered to have vanishing trilinear Yukawa 
couplings. This property does not follow from the gauge symmetries we considered and 
should be attributed either to discrete string symmetries or to additional gauge symmetries 
by enlarging the model. 6 

6. Neutrino physics 

One of the challenges of Standard Model extensions is the justification of the smallness 
of neutrino masses. The favorite scenario used to rely upon the introduction of right-handed 
neutrinos (SM singlets) and their mixing with some extra massive singlets. The suppression 
of the neutrino masses is then obtained as a result of the structure of the full mass matrix 
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("see-saw" mechanism). In order for this mechanism to work effectively, the extra singlet 
mass should be about ten orders of magnitude higher than the electroweak scale. 

Among the promising features of D-brane models is a novel scenario to account for 
neutrino masses: right-handed neutrinos are assumed to propagate in the bulk while left- 
handed neutrinos, being a part of the lepton doublet, live on the brane. As a result, the 
Dirac neutrino mass is naturally suppressed by the bulk volume. Adjusting this volume, so 
that the string scale lies in the TeV range, leads to tiny neutrino masses compatible with 
current experimental data. 

The extra dimensional neutrino mass suppression mechanism described above can be 
destabilized by the presence of a large Majorana neutrino mass term. As already mentioned 
in section 3, the lepton-number violating dimension five effective operator LLHH leads, 
in the case of TeV string scale models, to a Majorana mass term of the order of a few 
GeV. Even if we manage to eliminate this operator in some particular model, higher order 
operators would also give unacceptably large contributions, as we focus on models in which 
the ratio between the Higgs vacuum expectation value (VEV) and the string scale is just of 
order 0(1/10). The best way to protect tiny neutrino masses from such contributions is to 
impose lepton number conservation. As we have seen in section 2, we can find models which 
successfully accommodate all SM particles and preserve lepton number as an effective global 
symmetry in perturbation theory. These are the models A, A' and B, B' described in detail 
in section 3. 

Apart from neutrino masses these theories contain also the ingredients to explain neutrino 
oscillations. The right-handed neutrino, being a bulk state, has a tower of Kaluza-Klein 
(KK) excitations. Their mixing with the ordinary (left-handed) neutrino leads to oscillation 
patterns that have to be compared with present solar and atmospheric neutrino data. There 
exist extended discussions in the literature [3] regarding the neutrino mass and oscillation 
problems in the context of extra dimensional theories. Among the common results of these 
works is that an explanation of the solar neutrino anomaly is possible provided the Small 
Mixing Angle (SMA) solution is acceptable. However, recent SNO results in conjunction 
with SuperKamiokande data [221 Ell I2H I2H] strongly disfavor the SMA solution and thus 
render this higher dimensional oscillation mechanism problematic, at least as far as solar 
neutrino oscillations are concerned. A possible way out is to introduce three bulk neutrinos 
and explain the oscillations in the traditional way [2H] . The effect of the KK mixing can be 
eliminated by appropriately decreasing the size of the extra dimensions and thus increasing 
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the value of the string scale. However, all these discussions are restricted to the case of 
effectively one-dimensional bulk. Besides these phenomenological difficulties, there is also a 
serious theoretical problem, since one-dimensional propagation of massless bulk states gives 
rise to linearly growing fluctuations which yield in general large corrections to all couplings 
of the effective field theory, destabilizing the hierarchy j^J. 

Two-dimensional scenarios have not been considered in detail. We will see below how 
the above problems can be resolved and discover that a two-dimensional bulk has enough 
structure to describe both the solar and atmospheric neutrino oscillations by introducing a 
single bulk neutrino pair. On the other hand, recent experiments are also able to differentiate 
between the contributions of active and sterile neutrinos to the neutrino anomaly problems. 
From this point of view, the KK excitations do not carry any Standard Model charges and 
are thus considered as sterile. It is then important to examine if all these constraints are 
compatible with our model. 

As explained in the introduction, our setup incorporates a two-dimensional bulk and we 
are going to assume hereby that neutrinos propagate in the full bulk volume which is a 
two-dimensional space. Among the common features of the models considered in section 3, 
one finds tree-level neutrino couplings and mass-terms of the form: 

3 3 

^XiUHiVR -> ^2\iViU iL u R , (6.1) 

1=1 i=l 

where i is a generation index and for each generation i, Hi is one of the available Higgs 
doublets Hd or H u , providing masses to down quarks (models A, A') or to up quarks (models 
B,B'), respectively, with t>, = (Hi) the corresponding VEV. The above couplings provide a 
mass to one linear combination (i>l) of the weak eigenstates (va), while the other two remain 
massless. Note, that it would be possible to generate masses for all left-handed neutrinos 
by introducing additional bulk neutrino pairs. In this case the number of free parameters 
is increased and predictability is lost. Thus, here, we will study the case of a single bulk 
neutrino pair. Defining N L = [y L , v 0L , i/' 0L ) the mass eigenstates, the weak eigenstates can 
be written as 

ViL = U V Nl 3 ' ( 6 - 2 ) 

3 

where U is a 3 x 3 unitary matrix with = (i7 -1 )ij = and m? D = Y^h=i ^i v i * s the 
mass-square of the massive combination (vl)- Being of brane-bulk type, the couplings \ are 
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naturally suppressed by the bulk volume vgg (see section 1) and lead to a tiny Dirac neutrino 
mass 

v 2V2M S _ 

m D = —p== = -1T7-V , (6.3) 

Jv^ g 3 g 2 M P 



where v = y J^i=i hfvf with hi, i = 1,2,3, the associated dimensionless Yukawa couplings 
and Vi the corresponding Higgs VEV (v = (Hd), (H u )) depending on the model. Using 
typical values for the gauge couplings (see Table El of section 3), Wj < v — 246 GeV and 
hi/4-K = O (1), we obtain < 6 x 10~ 3 eV for M s ^ 10 TeV. This provides an explanation 
for the smallness of neutrino masses and is actually the extra-dimensional version of the 
see-saw mechanism. 

The above picture is simplified because we have neglected the contributions of the tower 
of KK neutrino states. Taking them into account, and assuming for simplicity that the two 
bulk radii are equal R$ = Rg = R and form an angle § — 0, where —tt/2 < 6 < tt/2, the 
mass terms become 



ml 
k 



L m = m D v h 4 fc) + E m k »R k) V R + c - c - ( 6 - 4 ) 

k k 

where Vr = v^) and the summation over k extends over all KK momenta. By m| we denote 
the mass-square of the KK excitation labelled by momenta k — (k\, k 2 ) 



ni 



2 



1 



k R 2 cos 2 e 



(k 2 + kl-2kxk 2 sm9) . (6.5) 



We also use the notation m\ for the mass-square of the fc-th KK level. 5 is a model dependent 
constant bigger than one, associated to the coupling of two Neumann-Dirichlet (ND) Z 2 - 
twisted strings to an untwisted (NN or DD) string [2Z1- I n our models, there are four 
Z 2 -twisted coordinates, implying 5 = 16. M plays the role of an ultraviolet (UV) cutoff, 
which is normally the string scale M s , but we prefer using the symbol M because in certain 
processes there exists an induced cutoff that can be a few orders of magnitude below M s . 
For instance, this is the case of solar neutrinos, where the production energy is of order of a 
few MeV, and thus heavier KK modes are effectively cut off. 

The mass terms (|6.4j) lead to a mixture of the usual left-handed neutrino with the infinite 
tower of its KK excitations. A detailed analysis of the eigenstate problem in our framework 
is presented in Appendix A where we derive the basic formulas for neutrino masses and 
transition probabilities. Due to its complexity, the problem can be either treated numerically 
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in the general case or analytically using some approximation. The first approach has the 
disadvantage of being rather tedious as it involves summations over a very large number 
of KK modes, so we will adopt here an analytic perturbative approach. Concerning the 
interpretation of neutrino anomalies, there are also two possible treatments: The first is a 
direct fit of neutrino data to the transition probability formulas obtained in our framework. 
The second is to try to simulate the standard solutions to the solar and atmospheric neutrino 
anomaly problems. We will use here the second method, as it is sufficient for demonstrating 
the basic features of our model. 

Following Appendix A, the mass spectrum of the full system, in the case of two bulk 
dimensions, is: 



"A 



ml = m z k + r k m 2 D 5 (1 - A fc ) + . . . (6.6) 

where r k is the multiplicity of the fe-th KK level and 

A k = irm 2 D (R 2 cos6) {log (M 2 R 2 cos 9 log 5 2 ) + s k } , (6.7) 

with s k a volume independent constant. Our solution is based on the assumptions that 
moR < 1, as justified by ()6.3|) and A n < 1 that simplifies the formulas involved. Under 
these assumptions, and following the analysis in Appendix A, the survival probability for a 
neutrino of flavor i is given by 

2/1 „ 2\ • 2f m oL\ 2 A „:„2 ( U L\ 



P v ^ Vi « 1 - 4<(1 - uf) sin 2 [~J - 3-2<A sin 2 f — — 1 , (6.8) 

where m 2 , = m 2 D {l — A ) and Ui = \Un\ satisfying the unitarity relation — 1; L is the 

distance that the neutrino travels before being detected and E is the beam energy. According 
to the discussion in Appendix A, the survival probability takes this form only for specific 
values of the angle 6: 

sm0 = ^, P ,qeZ, \p\<q ■ ^= \ t2 i , ( 6 - 9 ) 
q 2(g^ — p z ) 

where p, q are relatively prime integers and q = 1 for p = 0. In our approximation, the 
survival probability (|6.8j) is a superposition of two modes with frequencies: and j^^- 
These two frequencies can be considered as independent parameters, as the first depends on 
the Yukawa couplings and the Higgs VEVs, while the second depends on the compactification 
radius or equivalently on the string scale. The existence of these two frequencies provides us 



28 



with the opportunity to fit both solar and atmospheric oscillations using (jfTHJ). Furthermore, 
the amplitudes of the two modes depend on U{ and Ao defined in (16. 7|) . These parameters 
can be used in order to fit the oscillation amplitudes. 

In the standard neutrino (two flavor) scenario, one usually explains the solar neutrino 
anomaly by v e — > oscillations and the atmospheric neutrino deficit by — > v T oscillations. 
The formula for the transition probability is: 

'Am?- L 
~4 ~E 



sin 2 2% sin 2 ( — ) , (6.10) 



where Am?- = mf — m 2 is the neutrino mass difference in the case of two states mixing. 



Expressing L in kms, E in GeV and Am 2 in eV 2 , the frequency Am- takes the form 



1.27 x Am, 2 - L/E. 

Recent analysis of atmospheric neutrino data [24] at 3a CL. gives 1 x 1CT 3 < Am 2 tm < 
6 x 10 _3 eV 2 and 0.7 < sin 2 2# atm < 1. Regarding solar data, the situation has dramatically 
changed after the latest SNO results: Only the LMA and LOW MSW solutions are acceptable 
at the 3a CL. with 2.3 x 10~ 5 < Am 2 MA < 3.7 x 10~ 4 eV 2 , 0.6 < sin 2 2# LM A < 1 and 
3.5 x 10~ 8 < Aml ow < 1.2 x 10" 7 eV 2 , 0.8 < sin 2 2# LO w < 1. Moreover, the LMA gives a 
much better fit. The region of the SMA solution (with best fit values Amg MA ~ 5 x 10~ 6 eV 2 , 
sin 2 2#sma ~ 2x 10 -2 ) is acceptable only at the 5.5a level and is thus practically excluded [2*3] . 

The atmospheric neutrino oscillation frequency is higher than all solar ones, Am 2 tm > 
Am 2 ol , and thus we have to use the lowest frequency in ()6.8|) (i.e ml) to simulate solar neu- 
trino oscillations. Formula ()6.8|) contains four independent parameters, namely mc, R, M s 
and u e (assuming u T = and thus u 2 e + w 2 = 1). Fitting both solar and atmospheric 
oscillations requires to leading order in A : 



A : 






UJ 


= Am ltm 


(6.11) 




= Am sol 


(6.12) 


-<) 


= sin 2 26> so i 


(6.13) 


4^o 


= sin 2 20 atm 


(6.14) 



4u 2 e (l 



Neglecting the constant term so in the expression (j6.7j) of Ao, in the limit MR 3> 1, and 
assuming 5 = 16, Ao can be written in terms of M s and as 

1 m\Ml ( M P /bgT\ 2V2 
zna 2 M£ \ iraMg \ 2 I g 3 g 2 




where we have assumed that the cutoff is equal to the string scale (M = M s ). This choice 
of the cutoff is suitable for the atmospheric neutrino data, where the oscillation amplitude 
is proportional to Ao- In any case, the exact value of the cutoff plays a minor role in 
our calculation, due to the fact that it appears always logarithmically. Furthermore, the 
expectation value v is related to the rest of the parameters through equation ([6.3)1 . while the 
angle 9 enters in the Plank mass definition ()1.3|) 

1 M 2 

R 2 cos9= - = (6.16) 

In terms of the integers p, q that have been introduced in ()6.9|) . we can recastrewritw the 
last equation as 

W - W>' (6 . 17 ) 



or equivalently 



3 + (-l) q 

cos6 = — £A — L . (6.18) 



Thus, the four conditions (fB~TTll-(f6~Tll together with lj?TT5|l and (fOj) . (l6~TH|l fix all four 
parameters of the model. Therefore, fitting the atmospheric neutrino frequency ()6.11|) . one 
determines the compactification radius 



1 x l(T 3 eV 2 < -%r < 6 x 10~ 3 eV 2 , (6.19) 
R l 

or 3 /mi < R < 6 /jm for u; ~ 1. Choosing for the solar neutrino deficit the preferred LMA 
solution, we get from the second condition (|6.12j) the neutrino mass range: 

4.8 x l(T 3 eV <m D < 7.7 x l(T 2 eV . (6.20) 

The third condition (|6.13|) fixes the mixing coefficient u 2 e and has two possible solutions, 
namely, 0.18 < u\ < 0.5 or 0.5 < u 2 e < 0.82. Choosing u 2 e ~ 0.18 and u 2 ~ 0.82 (u T = 0), 
equation ()6.14|) leads to Ao ~ 0.27 (in the case we choose the lowest allowed value of 
sin 2 26 l a t m ), which lies at the edge of the validity of our perturbative approach. Any other 
choice of Ui compatible with the constraints leads to bigger values for A . This justifies also 
the choice u T = in order to minimize A in (|6.14p . 7 From (|6.14|) we get the string scale: 

8TeV£ M s £ 13TeV, (6.21) 



7 Normally, one should repeat the eigenstate analysis of Appendix A numerically in the non-perturbative 
region, but from a preliminary analysis we do not expect significant change of our results. 
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while compatibility with ()6.3|) requires 0(1) values for the Yukawa couplings. It is interesting 
that this range for the string scale coincides with the values we found from the analysis of 
gauge couplings in section 3, for the models B and B' . Coming to the angle, we get from 
(|6.17|) 0.02 < W p )(J < 0.2 for the allowed range of Am^ tm and we can easily verify that there 
exist integers p, q that satisfy (j6.17|) . 

Let us now consider the LOW solution to the solar neutrino deficit. Following similar 
steps, the four constraints (j6.11|) - (j6.14|) in this case give 

1.9 x 10" 4 eV <m D < 3.5 x 10~ 3 eV , (6.22) 

with u 2 e ~ 0.28 and thus u 2 ^ ~ 0.72 and Ao ~ 0.30 . The string scale turns out to be slightly 
lower in this case: 

1.8 TeV < M s < 2.2 TeV , (6.23) 

while for the angle 9 we get 20 < W p ^ q < 200. Note that the range of string scale is now 
compatible with the values found from the analysis of gauge couplings in models A and A' . 
In this solution, the left-handed neutrino Yukawa couplings do not have to be of O (1). 

Moreover, the practically excluded SMA solution can also be obtained in this framework. 
The associated parameters in this case are: u 2 e ~ 5 x 10~ 3 , u 2 — 1 — u 2 ~ 1, A ~ 0.2, 
m D ~ 2 x 10" 3 e1/, M s ~ 6 TeV, 0.2 < W PA < 1.2. Note that the case 9 = corresponds to 
p = 0, q = 1 and thus W VA = 1. As seen from our results, only the SMA solution includes 
this value in the allowed W-r&nge and this is the reason why only this solution could be 
reproduced in the case of an orthogonal torus. The LMA and LOW solutions require a bulk 
forming a non-orthogonal lattice, corresponding to non-trivial values of 9. It is also worth 
noticing that such non-trivial values of 9 induce CP violation in the neutrino sector, which 
is interesting to be further explored. 

The mixing of the neutrino zero mode with its KK excitations can lead to a decay of the 
left-handed neutrino to these KK modes, considered as sterile from the SM point of view. 
In our framework, and to leading order in the Ao expansion, the average conversion rate of 
a neutrino of flavor i to sterile is given by (jA.19|) : 

P Vi ^ a ~ 2u 2 A Q . (6.24) 

Constraints ()6.12|) and (|6.14j) fix both the above probabilities. Assuming the LMA solution 
to the solar neutrino deficit, we get 

P u ^ s ~ 0.44 (6.25) 
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for atmospheric and 

P v ^ s ~ 0.05 (6.26) 

for solar neutrinos, where in the second case we have assumed a cutoff M ~ 50 MeV. For the 
LOW solution, the transition probabilities are similar: P u ^ s ~ 0.32, P v<t -+ a ~ 0.08. Note 
that the decay rate to sterile neutrinos is significant in the case of atmospheric neutrinos and 
is negligible in the case of solar neutrinos. This is related to the structure of our model for 
neutrino oscillations. The atmospheric neutrino deficit is simulated using the lightest KK 
neutrino excitation (which is interpreted from the SM point of view as a sterile neutrino), 
while the solar data are explained using the (active form the SM point of view) zero mode. 

Constraints for the conversion of active to sterile neutrinos have been recently examined 
in reference [26J. Following their analysis in the case of the LMA solution, the constraint to 
the average decay rates for solar neutrinos is P Ve -+ Va < 0.40 at 90% c.l. which is obviously 
satisfied by our model. For atmospheric neutrinos the relative constraint takes the form 

AP = P u ^ Ua -P u ^ Us <0.17 . (6.27) 

Evaluating this constraint in our framework, one finds AP = 0.44 — 0.10 = 0.34 (where 
Pu e ^u s = 0.10 in the case of atmospheric due to the higher cutoff in A of eq (|6.7jl ) which 
is by a factor of two higher than the experimental bound. However, one should take into 
consideration that our perturbative analysis, focusing on explicitly revealing oscillations, 
does not allow to access the region A ~ 1 where in principle the above rates could change. 
As mentioned earlier this region could be studied only numerically. This requires summation 
over a huge number of KK modes and at present it appears insoluble even numerically. 

In any case the exact nature of atmospheric neutrino oscillations is expected to be further 
examined in the K2K [2E] experiment. In case the predictive scenario of a single bulk neutrino 
presented here fails to satisfy the sterile production constraints, one should proceed in the 
introduction of additional bulk neutrinos and explain oscillations in the traditional way, that 
is by zero mode mass difference and not by mixing with the KKs. Their presence can still 
lead to sterile production which can be reduced by appropriately raising the string scale and 
thus decoupling the KKs 
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7. Summary and conclusions 



In conclusion, we performed a systematic study of the Standard Model embedding in 
type I string theory at the TeV scale. We found that the minimum configuration with 
interesting phenomenological features requires three sets of D-branes, so that all SM particles 
are obtained as open strings stretched among these brane stacks. Two of them describe 
respectively the strong and weak interactions, while the third one contains a single abelian 
brane that extends in a two-dimensional bulk of submillimeter size. 

The model predicts the correct value of the weak angle for a string scale of a few TeV. It 
also contains baryon and lepton number as perturbative global symmetries, ensuring proton 
stability and absence of large (Majorana) neutrino masses. On the other hand, it uses two 
Higgs doublets that can provide masses to all quarks and leptons. Concentrating on the 
heaviest generation, we computed all trilinear Yukawa couplings and studied the resulting 
mass relations. We found a naturally heavy top and the mass ratio of bottom quark to tau 
lepton close to its experimental value. 

Finally, we have studied neutrino masses and oscillations by introducing a single right- 
handed neutrino state in the bulk. We found that both solar and atmospheric neutrino data 
can be explained if the bulk is a non orthogonal torus forming a non-trivial angle. Solar 
oscillations are then explained using the zero-mode, which obtains a tiny mass from the 
electroweak Higgs, while atmospheric oscillations use its first KK excitation. However, in 
the cases of atmospheric data, it seems to be an excess in sterile production with respect to 
current atmospheric data analyses. 

Overall, the model looks very promising and deserves further investigation. Particular 
directions that have not been discussed are the masses and mixing angles of the two lightest 
generations, possible important threshold corrections related to the two-dimensional bulk, 
supersymmetry breaking effects in models with brane supersymmetry breaking, as well as 
explicit type I string realizations. 
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Appendix A: Neutrino masses and oscillations 

We consider the neutrino mass eigenvalue problem that arises when the usual left-handed 
neutrino localized on the "brane" mixes with one pair of right-handed neutrinos propagating 
in a two-dimensional bulk. The solution of this problem in the case of one-dimensional bulk 
has already been studied in the literature j5]. A new feature of the two-dimensional bulk 
is that the associated KK sums are divergent and a mass scale M playing the role of the 
UV cutoff, normally identified with M s , appears in the mass and eigenstate expressions. 
Moreover, we consider neutrino oscillations and derive formulas for the transition rates of 
both active and sterile neutrinos. 

We will assume for simplicity that the two bulk radii are equal, Rg = R9 = R, but we will 
allow for the possibility that the angle 9 between the two compactified directions is arbitrary 
— 7r/2 < 9 < 7r/2. The masses of the KK excitations, labelled by momenta n = (jix, 712), are: 



— — (n? + n\ - 2n x n 2 sin 9) . (A. 1) 

R z cos z 9 v ' 

The KK modes can be ordered according to their mass and labelled by a unique level number 
k. Massive levels have in general degeneracy four, apart from particular points that have 
higher degeneracy for special values of 9. In any case, only the direct sum of the states of 
each degenerate level couples to the left-handed neutrino. Hence, we can diagonalize in the 
degenerate subspace and choose one of the eigenstates, which corresponds to the sum of the 
degenerate KK modes. In this basis, the relevant neutrino mass terms take the form 

2 

L m = m D v L ^2 \^kS~^ T + m k "r^ v^r + c.c. + decoupled , 

k k 

where = -^J2e u< R and ^r^ = 7^ Hi^r 1 ^ = 1> • • • > r k, with r k the multiplicity of 
the KK level with mass m^. The mass terms can be written in matrix form {Nj^mNn +c.c.) 
withiV L = (u L ,i>f\...),N R = (4 0) ,^ {1) ,...) andman infinite-dimensional matrix. 
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In order to determine the left-handed neutrino mass eigenstates, we consider 



mm) = 



/2 2 2 \ 

m 2 D A mxm D ^fr{b 'm 1 m^mo^fr^b m 7 ... m k m Dy /r^S m 7 » 



m\mr>^fr{b m 7 m\ ... 



m^mp^/r^b m 7 m| 



^/r k ~~im k _ 1 m D b 

m fc 

\ ^/TkmurnDb'^ 1 



ml 



(A.2) 



where A = X^ r ^~ ■ I 11 the sec L ue l we will assume that all masses are measured in string 
units, which we restore only at the end of our calculations. 

The exact eigenvalue equation for the mass of the n-th KK level m n can be written in 
the form 



and the associated eigenstates 

N, 



-2m? 



m^ — m. 



= «§,...) (A.4) 



with 



and 



c?=^^5-™2 (A.5) 



N 2 = m 2 D m 2 n V 1- . (A.6) 

u n ^ (ml - m 2 ^ 2 



The above results can be used to express v L in the basis of the mass eigenstates 



' N, 



n 



and calculate its time evolution 
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where E is the neutrino beam energy and L is the distance from the source. Therefore, using 
(|6.2|) we can derive the time evolution of the weak eigenstates 

va,(t) = UaVLit) + U l2 v 0L + U l3 v' QL . (A.9) 

The transition rate P Vi _> Vj , that gives the probability for a neutrino of a specific flavor i 
produced in the source to be detected as flavor j in the detector, is 



p,^ = \(^ L (o)WjL(t))\ 2 = < ' "v 1 . " (a.io) 



% = j 



a 



El / imiL\ 



where Wj = and 

T = Sr N* 2E 
The formulas for the transition probabilities to active neutrinos are 

p [ (1 - «?) 2 + « 2 (1 - «?) (T + T*) +ui\T\ 2 ,i=j 

{ ulu2(l-(T + T*) + \T\ 2 ) , i^j. 

Therefore, the transition rate for a neutrino of flavor % to decay into a sterile neutrino is: 

3 

P^ s = 1 - J2 P ^ = «?(!" \T\ 2 ) ■ (A.13) 

3=1 



Using (|A.11|) we obtain 



T + T* 2 ~mlL 

s i n 2 2_ o t? i^i 



1 - T72 sin -TET ~ 2F 2 (mi) (A.14) 
o 



T| 2 = l-4F 2 (ml)+4Fi(ml-ml) + F?(rhl) (A.15) 



where 



^) = E^ SinP #- ( A - 16 ) 

We can now calculate the average probabilities for a neutrino of flavor i to survive 

^ = (l-^ 2 ) 2 + « 2 C 2 , (A.17) 

or to be converted to flavor j 

P^=u 2 u 2 (l + ( 2 ), i^j, (A.18) 
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or to decay into sterile 

P v ^ s = ul{l-C 2 ) , (A.19) 

where 

n « 

and we have averaged over all frequency modes. 

For vtly) = 0, the eigenvalues of the matrix mm) form the usual KK tower with masses 
ma given by eq. (|A.1|) . For mjjR <C 1, it is natural to assume that the KK levels are slightly 
shifted: 

ml = m 2 n + 5m 2 n . (A.21) 

Inserting (|A.21|) into (|A.3|) and expanding for 5ml <C ml — mj V n ^ £, we obtain to lowest 
order: 

1 + A n ' { J 

where 

£-2mf <5~ 2m f 

Due to the presence of the factor <5~ 2m «/ M2 (after restoring the M units), we have 5ml ~ 
for m n > M, implying that KK levels above the cutoff M are not shifted. We notice also 
that in (jA.lj) we have R2 * os g ~ -^f- <C M 2 . In this limit, we can calculate the leading terms 
in A n : 

A n = vr m 2 D (R 2 C os 9) {log (M 2 R 2 cos 9 log 5 2 ) + s n } , (A.24) 

where s n is a constant term independent of the cutoff. Specifically, Sq = C and s n -^o = 
— log \n\ + n 2 — 2riiH2 sin9\ + C, with C a constant of order one. Similarly, for the normal- 
ization coefficients we get to the lowest order in 5m 2 n 

1 Sm h... (A.25) 



N 2 m D 



N 2 m 2 D m 2 n r n 5 2m " 
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(5m « )2 + . . . (A.26) 



The infinite KK sums F2, F\ in (|A.14|) . ()A.15|) can in principle be calculated numerically 
using (|A.3|) and (|A.6|) . However, this requires summation over a huge number of KK modes. 
A convenient approximation is to assume A„ << 1 for n > 0. In this region 

n^(0,0) 

and thus, F 2 can be expressed in terms of theta-functions using the formula 

f - 2m 2sin 2 {13ml) f r / 4(xlog5 + i) \ ( 4(xlog5 + z; 



^ rr> 2 L 1 6 \ tt( -I — ai-n ftM?2 



ml J I °V tt(1 - sin 6) R 2 J °\ tt(1 + sin#)i? 2 

where $3(1") = Xlnez e 271 "™ 2 and #2 ("7") = Snez e w v +5 ) . In the case of an orthogonal torus 
(9 — 0), F is periodic under ^ — > ^ + 7ri? 2 . However, this property is lost for arbitrary 
values of 9. The periodicity, of the survival probability (|A.10jl . is necessary for interpreting 
the neutrino anomaly through neutrino oscillations and is in general restored for rational 
values of the angle 9 

sin0 = - , \p\<q, (A.29) 
Q 

where p, q are relatively prime integers, and q = 1 for p = 0. Restricting 9 to this subspace, 
F becomes periodic under — > + T P!q iiR 2 , where 

{- — for q = odd 
f-v> , < A ' 30 > 
for o = even. 

Since F 2 is a periodic function, the next question is to compute the corresponding am- 
plitude. To get an estimation of the amplitude we can evaluate the sum (jA.27|) at the half 
period. We get 

r— 2m- 

^ 1/2 = ^dE^ = kA 0' ( A - 31 ) 
z — ' mt 

where Z 2 is one of the sets of (even, odd) and/or (odd, odd) integers of Z 2 , depending on 
the choice of q. More particularly, for q = 4£ with £ integer, Z 2 is the set of (odd, odd) 
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pairs, for q = 2£ + 1, Z 2 is the set of (odd, even) and (even, odd) pairs and for q = 4£ + 2, 
Z 2 is the union of the two previous sets. The constant k takes approximately the values 
k G (1/4,1/2,3/4) for each of the three cases, respectively. Moreover, an upper bound to 
the amplitude of F can be derived by replacing the "sin 2 " terms with unity 

F 2 max = m 2 D V — ^ = A . (A.32) 

1/2 

Hence, the oscillation amplitude p lies in the range F 2 < p < F™ ax , that is kA < p < A 
with k G (1/4,1/2,3/4). Furthermore, we can proceed to a numerical evaluation of F for 
given sin 9 and MR. An explicit example is presented in Figure where we have calculated 
F as a function of in the case M 2 R 2 = 10 4 , 5 = 16, sin^ = In the same figure, we 
have also plotted the function sin 2 (H§ 4E#0 (S Ta J nne ) with amplitude arising from the fit 
to the numerical sum data. The numerical evaluation of the sum shows that the amplitude 
p can in general be approximated by p ~ 0.8 Ao. Taking into account these results, we will 
assume in section 6 that the KK sum F can be simulated by the dominant frequency mode 

F,«0.8W(^) ; ^f^. (A.33) 

where Ao is given in f)A.24|) and 6 is given in ()A.29j) . 

Assuming A 1, we can drop the terms F 2 2 and F 2 < F 2 2 in (jA.15|) . Putting together 
(jA.12|) . (jA.27|) and (jA.33|) . we obtain an approximate expression for the survival probability: 

« 1 - 4 Mi 2 (l - uf)(l - A ) sin 2 ^ - 3.2 M 2 A sin 2 , (A.34) 

where 

m 2 = m^(l-A ). (A.35) 

Moreover, for the parameter ( 2 that enters in the average transition probability formulas 
(|A~T7)l . (lA~T8l) . (|A~T9)l . we have 

C 2 ~ = 1 - 2A + . . . . (A.36) 
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